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Abstract 

We present an atomistic investigation of the influence of strain on the electronic properties of 
quantum dots (QD's) within the empirical sp 3 s* tight-binding (ETB) model with interactions 
up to 2nd nearest neighbors and spin-orbit coupling. Results for the model system of capped 
pyramid-shaped InAs QD's in GaAs, with supercells containing ~ 10 5 atoms are presented and 
compared with previous empirical pseudopotential results. The good agreement shows that ETB is 
a reliable alternative for an atomistic treatment. The strain is incorporated through the atomistic 
valence force field model. The ETB treatment allows for the effects of bond length and bond 
angle deviations from the ideal InAs and GaAs zincblende structure to be selectively removed from 
the electronic-structure calculation, giving quantitative information on the importance of strain 
effects on the bound state energies and on the physical origin of the spatial elongation of the wave 
functions. Effects of dot-dot coupling have also been examined to determine the relative weight of 
both strain field and wave function overlap. 

PACS numbers: 73.22.-f, 68.65.Hb, 71.15.Ap 
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I. INTRODUCTION 



Nanometer-size semiconductor quantum dots (QD's) have attracted scientific interest due 
to their potential applications in optoelectronic devices as well as because of their peculiar 
properties such as the self-assembly, tunability, narrow size distribution and large Coulomb 
blockade effects.— The size and shape of the Stranski-Krastanow growth of InAs QD's on 
GaAs (001) reported by different authors vary, depending on the epitaxial method and on 
the growth conditions. Different sizes of QD's, pyramidal or dome shapes with side facets 
oriented along different directions,— i 3 i 4 i 5 i 6 truncated cone£ and pyramids^ with nonuniform 
Ga incorporation in the nominally InAs QD's have been reported. The driving force for the 
formation of such structures is the relief of the elastic energy associated with a dislocation- 
free, epitaxial structure (the InAs/GaAs lattice mismatch is 7 %). The strain distribution is 
not uniform, so accurate electronic models should include the effects of such nonuniformity. 

Theoretical models currently employed in the study of the electronic properties of QD's 
can be generally divided into macroscopic or microscopic. Examples of macroscopic models 
are the one-band effective mass approximatio n 9 ' 10 ' 1 1 and the multi-band k ■ p models . 12 ' 13 Mi- 
croscopic models are based on the empirical pseudopotential method^ and on the empirical 
tight-binding (ETB) method . 15 ' 16 ' 17 ' 18 ' 19 ' 20 ' 21 ' 22 ' 23 Microscopic models provide an atomistic 
treatment, as required for a more realistic description of smaller heterostructures. Here, 
the effects of inhomogeneous strain follow directly by taking into account deviations of the 
atomic positions from the ideal InAs and GaAs bulk structures. The empirical pseudopo- 
tential treatment potentially offers the most accurate description of the electronic properties 
of QD's. On the other hand, the ETB method may offer a faster alternative, and it is more 
transparent with respect to the analysis of results in term of the modified chemical bonding 
present in and around quantum dots. 

Despite their potential strength, not many ETB calculations are available regarding 
capped strained InAs QD's. The only published results are obtained by modeling the dot by 
a spherical cluster with dangling bonds saturated by hydroge n 18 ' 19 ' 20 ' 21 (this approximation 
is valid in the limit of nanocrystals embedded in a material with a very wide gap) or by a 
pyramidal dot with uncovered surfaces.— 

In this work we explore the ETB method for evaluating and analyzing the electronic 
structure of InAs QD's. Our aim here is to: 
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1. Study the reliability of the ETB scheme for the treatment of semiconductor nanos- 
tructures. We consider a square-based pyramidal InAs QD with {101} side facets 
embedded in GaAs. Since this geometry has been previously adopted by several 
authors , 12 i 13 i 14 i 17 the reliability of our results is assessed through comparison with 
previous studies. 

2. Investigate how the strain affects the electronic properties of QD's. So far, such in- 
vestigations have been limited to spherical 2 ^ or elliptical 2 ^ dots. In the first case, the 
influence of the strain was estimated by comparing free standing with GaAs embedded 
quantum dots. However the surface dangling bonds in the free-standing dot were pas- 
sivated by a fictitious material with a band gap much larger than the GaAs gap, giving 
raise to a much larger confining effect for electrons and holes inside the dot. Therefore 
in the comparison not only the different strain configurations played a role, but also 
the different band offsets. In the second case a comparison of results of differently 
strained dots was made, without however including any bond angle deformation. In 
the present study we complement these results by exploiting the flexibility of the ETB 
formalism where strain effects may be entirely removed from the model Hamiltonian 
(without any structural simplifying assumption), allowing direct comparison of the 
real QD with an artificially strain- unaffected GaAs-embedded QD. 

3. Analyze the influence of strain field and inter-dot hybridization for well separated dots. 
Previous studies of inter-dot couplin g 11 ! 26 focused on a complementary range of dot 
separations, namely closely stacked dots. 

We calculate the single-particle bound electron and hole states and wave functions adopt- 
ing Boykin's sp 3 s* parametrization with interactions up to 2nd nearest neighbors and spin- 
orbit coupling.— This parametrization gives very good fits of the important effective masses 
and gaps for bulk GaAs and InAs. One potential problem is that it does not reproduce 
as well the d-bands contributions as parametrizations that explicitly include ci-orbitals, as 
for instance the one proposed by Jancu et a/.,— where a sp 3 d 5 s* basis set and 1st nearest 
neighbor interactions were considered. Since the electron bound states in our QD come 
mainly from s and p atomic states, this does not constitute a relevant drawback. We do 
not consider piezoelectric effects in our model. However, as we will remark in Sec. IIII Bl 
we expect only minor corrections to our results coming from such effects. The influence of 
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strain on the bond lengths is taken into account through a power-law scaling of the ETB 
parameters chosen here as such as to reproduce hydrostatic pressure effects in both bulk ma- 
terials, while the influence of strain on the bond angles is taken into account by a generalized 
Slater-Koster formalism.— 

This paper is organized as follows. In Sec. Owe present the formalism adopted for struc- 
tural relaxation of the system as well as for the electronic structure calculations, including 
the geometrical power-law scaling of the ETB parameters. Our results are given in Sec. IHIl 
and in Sec. IIVI we present a summary and conclusions. 

II. FORMALISM 

A. Structural analysis 

In order to calculate the atomic structure of an InAs quantum dot embedded in a GaAs 
matrix, i.e. the strain relaxation, two different methods have been used in the literature. One 
approach is an extrapolation of continuum nonlinear elasticity (CE) theory to the atomic 
scale, employing a discretization which is either based on the finite differences^^ 3 - or the 
finite-element (FE) method . 30 ' 31 The alternative is the valence-force field (VFF) approach, in 
particular the Keating model . 32 ' 33 ' 34 The latter approach has several advantages : it accounts 
for internal displacements between the two sublattices of a zincblende crystal, that cannot 
be addressed within conventional continuum elasticity theory, and gives a displacement 
field which obeys the correct symmetry group Ci v - 35 However, for large systems and slowly 
varying strain fields, the computational effort using the VFF approach is higher than using 
the FE method because, in the FE calculation, atomic resolution is usually not required 
in all regions of space. In the present study, we started from the continuum elasticity 
theory as implemented in the FE method (using experimental elastic constants 36 ) to get 
a first approximation to the displacement fields. Then, by interpolation for the atomic 
positions that lie between the nodes of the FE's, we extracted all the positions of the 
atoms in the supercell that give rise to the calculated displacement field. Finally we refined 
these displacements by further relaxing the atomic positions using a VFF model. Thus the 
positions of the atoms are eventually determined within the VFF model, thereby ensuring 
the correct C 2v symmetry of the atomic displacements fields. 
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In the VFF model, the elastic energy energy of a zincblende lattice is expressed as a 
function of the atomic positions {Rj} as 



Here, <i°- denotes the bulk equilibrium bond length between nearest neighbor atoms % and 
j in the corresponding binary compound, and 9q = cos _1 (— 1/3) is the ideal bond angle. 
The first term is a sum over all atoms i and its four nearest neighbors j, the second term 
is a sum over all atoms i and its distinct pairs of nearest neighbors j and k. The local- 
environment-dependent coefficients, and (3^ are the bond-stretching and bond-bending 
force constants respectively. We use 3 ^ for GaAs = 2.448 A, = 41.49 x 10 3 dyne/cm, 
f3 ijk = 8.94 x 10 3 dyne/cm; for InAs we use d£. = 2.622 A, a = 35.18 x 10 3 dyne/cm, 
Pijk — 5.49 x 10 3 dyne/cm. Across the heterointerfaces, where the species j and k are 
different (Ga and In), we use for f3 the geometrical average of the corresponding values for 
pure GaAs and InAs. 

The elastic energy is minimized with respect to the atomic positions {Rj }. In the 
minimization process, each atom is moved along the direction of the force on it, Fj = — Vj-E 1 , 
and the movement is iterated until this force is smaller than 0.001 eV/A. 

We compared the elastic constants derived from the VFF model to the experimental 
ones. The elastic constants G\\ and Cyi agree with the experimental values within 6 %. 
Differences are noticeable mainly in C44. The VFF model gives the C44 about 10 % too 
low for GaAs and about 20 % too low for InAs. In order to estimate the error due to 
inaccurate elastic constants, we calculated the local strain tensor by CE using both the 
elastic constants derived from the VFF model and the experimental ones. By comparing 
these results we verified that the absolute error in the diagonal components of strain tensor 
was always smaller than 0.005. 
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B. Electronic calculations 



The electronic structure is obtained within the ETB approach, adopting a sp 3 s* 
parametrization with interaction up to 2nd nearest neighbors and spin-orbit coupling,— 
which has been successfully used for III-V semiconductor heterostructures . 37 ^ 38 The wave 
functions are written as \1/ = Ylwa c iva\iv ><j, where \iv > a are orthogonal normalized 
Ri-centered orbitals of angular type v = s,p x ,p y ,p z , s* and spin a, and c\ va are complex 
expansion coefficients. The s* orbital was first introduced by Vogl et a/.— to obtain a better 
description of the conduction bands. In a iV-atoms system, the lOiV x lOiV ETB Hamilto- 
nian matrix contains 33 independent matrix elements for bulk GaAs and 33 for bulk InAs. 
These matrix elements are the parameters of the model for the present calculation, and are 
taken from Ref. In a strained InAs/GaAs mixed material, such as the QD system, a 
new parameter related to the valence band offset also needs to be included in the model. 
This parameter consists in a shift of all diagonal Hamiltonian matrix elements for bulk InAs 
(resulting in an analogous shift of the InAs bands), and it has been chosen such that the 
energy difference between the bulk InAs and the bulk GaAs valence band edge to coincide 
with the bulk valence band offset (A„). 

We performed an analysis of the QD gap dependence on the specific choice of A„, because 
there is a considerable spread in the experimental values reported for A„ in the literature.— 
By varying A„ in the range 52 - 300 meV we obtained a QD gap variation smaller than 4 
%, indicating that in this range our results are not much affected by the specific choice of 
the offset. In what follows we take A v =52 meV from Ref. 3, in order to better compare 
with results reported there . 

The relaxed geometry of the QD system implies changes in bond lengths and in bond 
angles as compared to the ideal bulk materials. Both effects are incorporated in our electronic 
model. Bond length deviations with respect to the bulk equilibrium distances g??- introduce 
corrections to the ETB Hamiltonian off-diagonal elements 14/ ■ Note that recently a different 
scheme has been proposed,— where corrections to the diagonal matrix elements have also 
been included, in the framework of the sp 3 d 5 s* 1st nearest neighbors parametrization.— We 
assume a power-law scaling^ for the off-diagonal elements : 




(2) 
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TABLE I: Volume deformation potentials (in eV) for direct (a^ _r ) and indirect band gaps (a 
and a^~ x ), absolute volume deformation potential for the conduction band edge (a£ 6c ) and defor- 
mation potentials for uniaxial strains along [001] (b) and along [111] (d) (see text). Uniaxial strains 
were applied starting from the experimental 3 ^ lattice constants. For the uniaxial strain along [111], 
the internal atomic displacement is calculated using the VFF method. Our ETB calculations are 
compared with experimental results 36 (Exp.), with DFT-LDA calculations using pseudopotentials 
(PP)44 and with DFT-LDA calculations using the LAPW method (LAPW) & 



where |R, — Rj| is the actual bond length and Vki(dfj) is the bulk matrix element taken from 
Ref. |27| (k and / label the different matrix elements). The exponents Uki are determined to 
reproduce variations of the relevant binary materials electronic properties under hydrostatic 
pressure, namely the volume deformation potentials a" 



a a = V — ^ 
OV 



(3) 



V 



for a corresponding to the direct as well as indirect (T — L and T — X) band gaps. 

In Table H] we give the values for a" for GaAs and InAs taken from experiments— from 
a density-functional theory (DFT) calculation using the local-density approximation (LDA) 
and ab initio pseudopotentials (PP),— from a DFT-LDA calculation with the linearized 
augmented planewave (LAPW) method,— and from our results (ETB). In principle each nki 
depends on the orbital character. However, we find that a single exponent n^i = 3.40 for all 
integrals and both materials gives a satisfactory agreement with LAPW for a^~ r and a^~ L . 
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For a^~ x the agreement is less satisfactory. We believe that this difference reflects the fact 
that the bottom of the conduction band at X has a noticeable d contribution.— Therefore the 
inclusion of d states in the parametrization and a correspondent different Uki value (cf. Eq.|2J) 
would be necessary. However for InAs/GaAs QD's at atmospheric pressure, the confinement 
effect for electrons and holes inside the dot comes from the conduction and valence band 
offsets at the V point,— whereas the X point does not play an important role. Therefore 
we do not expect this disagreement to be relevant in our calculations. It is interesting to 
note that this single exponent is very close to the value n = 3.454 reported for GaAs and 
AlAs within a different ETB parametrization.— In Table H] we also give the absolute volume 
deformation potential for the conduction band edge (a^ 6c ), and the deformation potentials 
for uniaxial strains along [001] (b) and along [111] (d), obtained by^ 



6e 001 = 3 6(C-0 

5e in = 3V3 d elf (4) 

where <5eooi and 5em are the energies of the light hole band with respect to the heavy hole 
band, in the absence of spin-orbit coupling, for strains along [001] and [111] respectively, 
and e^ 1 and ejj 1 are components of the correspondent strain tensors, defined as 
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where Cn, Cxi and C44 are the experimental? 6 elastic constants. For the uniaxial strain 
along [111], the internal atomic displacement is calculated using the VFF method. 

Fig. H confirms the adequacy of the single-exponent scaling for the present study by 
comparing our ETB with DFT-LDA results for the InAs and GaAs T — T band gaps obtained 
by varying the lattice constant. The DFT-LDA calculations were performed using scalar- 
relativistic ab initio pseudopotentials of the Hamann type.— The electronic wave functions 
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FIG. 1: Comparison between ETB and DFT-LDA results for the InAs and GaAs r — T band gaps 
obtained by varying the lattice constant. We also report the T — X band gap calculated by ETB. 
The vertical dashed lines mark the bulk lattice constants (6.055 A for InAs and 5.653 A for GaAs). 

were expanded into a plane-wave basis set with a cut-off energy of 16 Ryd. It is well 
known that the band gap is underestimated in LDA, but the overall behavior of the gap vs. 
hydrostatic lattice deformation should be reliable. Our ETB scheme yields the experimental 
optical band gap and the volume deformation dependence follows closely the trend obtained 
with the LDA for a wide range of deformations. On the same figure, we also reported the 
r — X band gap calculated by ETB. We can observe that the band at the X point has 
higher energy than at the V point, for the whole range of lattice distortions typical in a 
QD (where the InAs (GaAs) is compressed (expanded) by at most 7 %). This behavior 
would not change if the ETB-calculated a^~ x reproduced better the LAPW results (see 
Table H}, since the curve representing the T — X band gap has a slope much smaller than the 
T — r band gap curve, and the crossing point between them would not change its position 
appreciably. These considerations confirm that a more accurate a^~ x would not affect the 
results presented here. 

Bond angle distortions are included in the ETB Hamiltonian as suggested in the Slater- 
Koster formalism,^ generalized to include three-center integrals for the 18 independent 2nd 
nearest-neighbor matrix elements.-^ Note that, different from previous studies^, we do not 
assume that the three-center integrals are independent of directional changes induced by the 
strain. 
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The relevant eigenstates of the resulting Hamiltonian matrix including all the strain 
effects, if, are obtained variationally. We build the quotient 



where e r is a reference energy. By minimizing 3? with respect to the trial function ip by a 
steepest descent algorithm, we get the eigenvector (and the related eigenvalue) whose energy 
is nearest to e r . Therefore by varying e r we may in principle determine all the electron and 
hole bound state energies and wave functions . 37 i 38 i 47 i 51 

Within our ETB formalism, strain effects can be formally removed from the electronic 
calculation by imposing = in Eq. (J2J) (removal of the strain and relaxation effects 
from bond lengths) and setting the direction cosines between atomic orbitals equal to the 
corresponding bulk values (removal of the strain and relaxation effects from bond angles). 
Therefore, by contrasting the bound state energies of an artificially strain-unaffected QD 
with the corresponding results for the strained QD, we are able to quantify the total strain 
impact on the electronic properties. 

III. RESULTS 

A. Relaxed QD geometry 

Fig. El is a schematic view of our pyramidal InAs QD buried in a GaAs matrix. The 
wetting layer is modeled by a monolayer-thick InAs layer at the base of the pyramid. The 
pyramid base length is 12a, the height is 6a, where a=5.653 A is the lattice constant of 
bulk zincblende GaAs. We place the InAs pyramid and wetting layer in a large GaAs box, 
to which periodic boundary conditions are applied. This supercell is used in our structural 
and electronic calculations. Three different supercells were considered, containing the same 
QD but differing by the size of the GaAs matrix, namely GaAs matrices with dimensions 
19a x 19a x 14.067a (40 432 atoms), 25a x 25a x 17.067a (85 000 atoms - shown in Fig. |2D, 
and 37a x 37a x 35.067a (383 320 atoms). The z-dimension is not an integer multiple of a 
due to the InAs wetting layer. Unless specified otherwise, results presented below refer to 
the 85 000 atoms supercell. 

In Fig. El we present the relative distortion of the bond angle 9 from the ideal zincblende 
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14.1 nm 




.6 nm 



FIG. 2: Schematic view of the pyramidal InAs QD buried in the GaAs matrix. The wetting layer 
is one monolayer- thick. The pyramid base length is about 6.8 nm, the height is about 3.4 nm. 
The supercell contains 85 000 atoms. 




FIG. 3: Average relative distortion of the bond angle from the ideal zincblende bond angle 6q, 
and average relative distortion of the bond length d from the ideal GaAs and InAs bond lengths do- 
These quantities are calculated along the [110] direction at z = 0.4 h. For comparison we present 
results obtained from the 40 432-atoms supercell (Small) and the 383 320-atoms supercell (Big). 

bond angle 6q (obtained by averaging over the 6 different bond angles around each cation), 
and the relative distortion of the bond length d from the ideal InAs (inside the dot) and GaAs 
(outside the dot) bond lengths do (averaging over the 4 bond lengths). These quantities are 
calculated along the [110] direction at z = 0.4 h, where h is the QD high. We compare 
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results obtained from the 40 432-atoms supercell (Small) and the 383 320-atoms supercell 
(Big). The crystallographic directions are defined by starting from a zincblende unit cell 
containing a cation at the origin and an anion at a(|,|,j)- The figure shows that the 
thicker GaAs region allows for a better relaxation of the bond lengths in the largest part of 
the supercell, while it is not as efficient in angular relaxations. 

B. Bound states of the relaxed QD 

We have calculated electron and hole bound states, and refer to them as |el > and \hl > 
for the respective ground states, |e2 > and \h2 > for the next excited states, and so on. For 
the QD supercell containing 85 000 atoms the energies are calculated as 1405 meV (|e2 > 
state), 1267 meV ( | el > state), 74 meV (\hl > state) and 39 meV (\h2 > state) using the 
top of the bulk GaAs valence band as energy zero. These energies are shown in Fig. EJ on 
the left side (column labeled QD). From our numerical approach, we cannot exclude the 
possible existence of other hole states with smaller energies and very close (Ae < 10 — 15 
meV) to the \h2 > state. 

We show in Fig. on the left side (QD), the isosurface plots of the charge densities 
e|yj(r)| 2 corresponding to the electron states |e2 > and | el >, and to the hole state \hl >. 
The isosurfaces are selected as 0.5 of the maximum charge density value. The figure shows 
that the charge is almost entirely confined inside the dot. The lowest electron state ( | el >) 
is almost s-like (slightly elongated along [TlO]), the next electron state (|e2 >) is p-like 
aligned along [TlO], and the hole state (\hl >) has an elongation perpendicular to the | el > 
state, in agreement with the work of Stier et alM and Wang et alM In Table |H] we show 
a comparison of the energy differences between the QD bound states calculated within the 
empirical pseudopotential approach (PP) (whose results are extracted from Fig. 2 of Ref.Q) 
and the present approach (ETB strained). We can see that the agreement is good. 

In larger pyramidal QD's, an additional p-like electron state oriented perpendicularly to 
|e2 > is usually present . 13 1 14 When calculated within macroscopic models, these two p-like 
states are degenerate, if the piezoelectric effects are neglected.— In our calculations we did 
not find this additional p-like state, and we observed that the |e2 > state lies about 20 
meV below the GaAs conduction band edge. On the other hand, when we considered the 
artificial strain-unaffected QD (see next section), all the electron states became deeper, and 
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QD QD QD 

strain-unaffected strain distances 

1 424 = GaAs CBM 

I— 1382— |e2> 

/ — 1214— |el> 
1155—' / 

948 — / 

FIG. 4: QD bound state energies (in meV) calculated by our ETB approach. The energy zero is 
the bulk GaAs valence band maximum. Different degrees of strain are taken into account. In the 
first column (QD) the bound state energies of the physical QD are reported, where the strain effects 
have been included in the ETB hamiltonian. The second column (QD strain-unaffected) gives these 
energies for an artificially strain-unaffected QD, discussed in Sec. IIIIGI The third column (QD 
strain distances) gives results when strain is retained only in the bond length description, while 
bond angles are assumed to equal the bulk ones (discussed in Sec. IIII G|) . 



PP ETB 

£|e2> - £|el> 130 138 

e lel> -e [hl> 1150 1193 

e \hl> — e \h2> 25 35 



TABLE II: Comparison of the QD bound state energy differences (in meV) obtained from empirical 
pseudopotentials^ (PP) (see text) and the present ETB results (ETB). 

the additional p-state did appear about 30 meV above |e2 >. This degeneracy lifting appears 
in our atomistic model as a consequence of the breaking of the pyramidal C± v symmetry into 
the lower zincblende C2 V symmetry. 

We now analyze how the first electron and the first hole states are affected by the super- 
cell size. Due to the periodic boundary conditions, different cell sizes correspond to periodic 
three-dimensional QD arrays with different inter-dot separations. In each case, before per- 
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FIG. 5: Isosurface plots of the charge densities e|(/?(r)| 2 relative to the electron states |e2 > and 
| el > , and to the hole state \hl >. Each surface correspond to 0.5 of the maximum charge density 
value. The left side (QD) shows the results obtained from the physical QD, while the right side 
(QD strain-unaffected) shows the results obtained from the strain-unaffected dot (discussed in 
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FIG. 6: Single-particle electron and hole energies (e| el> and €|/ ll> respectively) and QD gap (e| el> — 
e|/ ll> ) as functions of the distance between dots along [001]. The solid curves are phenomenological 
exponential fits, while the horizontal dashed lines are their asymptotes. 

forming the electronic calculation, the atomic positions are relaxed as described in Sec. Ill Al 
Supercell-size effects are due to electronic and elastic dot-dot interactions. Both contribute 
to the results shown in Fig- El where the energies for the states |el > and \hl > and the QD 
gap are shown for the three different supercells. The horizontal axis represents the supercell 
dimension along [001], i.e. the base-to-base inter-dot distance along the [001] direction. Al- 
though we have chosen to report our results here and in the following section as a function 
of the inter-dot distance along this direction, dot-dot interaction effects in all directions in 
the three-dimensional QD array are included. We note that the dot-dot coupling in the 85 
000-atoms supercell makes the gap wider by about 12 meV with respect to an isolated dot. 
Strictly speaking, when we bring QD's together to form a periodic array, the bound states of 
the isolated dot spread into minibands whose width increases with the dot-dot interaction. 
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Now a brief remark on the possible influence of the piezoelectric effects on the results 
shown in Fig. El (and in the subsequent Fig. [TJ, in the next section). Grundmann et a/.— have 
shown that the piezoelectric potential inside a pyramidal InAs/ GaAs QD's has a quadrupole- 
like character in the [001] plane. Moreover Fig. |5] shows that our | el > and \hl > states 
are almost symmetric under rotations around a [001] axis passing through the tip of the 
pyramid. It follows that in the framework of nondegenerate perturbation theory (where the 
piezoelectric potential is taken as perturbation), the first-order corrections to the energies 
e| e i> and e\hi> almost vanish. We therefore expect that the inclusion of the piezoelectric 
effects would not strongly affect the results shown in Figs. El and 

C. Bound states of an artificially strain-unaffected QD 

The influence of strain on the electronic properties was studied here by comparing the 
physical QD bound states with the corresponding artificially strain-unaffected QD states, 
as explained in Section Til Bl The second column in Fig. 0] (QD strain-unaffected) gives the 
bound state energies for the artificially strain- unaffected QD, that should be compared with 
results from the physical dot (QD). Note that the energy e\ h2> for the strain-unaffected dot 
is not given, because the state \H2 > is unbound. We observe that strain increases the QD 
gap (e[ei> —€\hi>) by about 25 %, raising it from the strain-unaffected value 937 meV to the 
value 1193 meV. This behavior comes mainly from the InAs main gap increase when the 
structure is compressed by the surrounding GaAs matrix,— as can be seen in Fig. ^ in the 
case of bulk InAs. The change in the QD gap due to strain reported here agrees qualitatively 
with effective mass calculations in elliptic dots.— We note that the band-gap variation with 
the strain (wider or narrower gap) is dominated by the |el > state position (shallower or 
deeper confinement). We observe that the strain has opposite effects on electron and hole 
states: electron states become shallower, approaching the conduction band edge, while hole 
states become deeper, moving far from the valence band edge. 

The last column of Fig. 0] (QD strain distances) gives results obtained by retaining in the 
Hamiltonian only the bond length deformations, assuming bulk bond angles. This shows 
that the electronic properties of the QD are mainly affected by deviations of the bond lengths 
from the respective bulk ones. We note that the QD gap for the "partially strained" system 
(1211 meV) is larger than for the physical QD. In fact the hole level actually drops with the 
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bond length compression, but it rises with the bond angle distortions resulting from the QD 
geometry. The angular contribution dominates, leading to the smaller gap in the physical 
QD case. We also observe that the angular strain contribution is more important for \hl > 
(71 meV) than for \el > (53 meV), since the former has wave function atomic components 
predominantly p-type, while the latter has wave function atomic components almost purely 
s-like (thus spherically symmetric). 

The right hand side of Fig. E] (QD strain-unaffected) shows the isosurface plots of the 
charge densities of the strain-unaffected QD. By comparing with the results of the physical 
dot, we observe that the spatial orientation of |e2 > does not depend on the mesoscopic 
Civ symmetry (resulting from the strain field), but depends on the alternating interface 
structures of the four {101} facets (resulting from the microscopic zincblende structure). On 
the other hand, the spatial elongations of |el > and \h\ > depend only on the mesoscopic 
Civ symmetry of the strain field, and not on the alternating interface structures of the four 
{101} facets. 

In Sec. IIIIBI we discussed the supercell-size effects on the electronic and elastic dot- 
dot interactions. The results shown in Fig. El reflect the effects of both interactions. By 
repeating now the calculation for the | el > and \h\ > states in the three strain-unaffected 
supercells, we were able to isolate the electronic interaction due to the wave functions overlap. 
The results are shown in Fig. [TJ where the solid curves are phenomenological exponential 
fits and the horizontal dashed lines their asymptotes. All the fits are of the form e = 
eo + A ■ exp(—d/\), where d is chosen as the base-to-base inter-dot distance along [001], eo 
represents the energy of the isolated dot (given by the horizontal line), A is the characteristic 
length of the interaction along [001], and A is a prefactor related to the interaction. For 
the overlap contribution (Fig. [TJ), an exponential dependence is to be expected since this is 
the typical behavior of the localized wave functions away from the dot. For the strain field 
contribution, a power law dependence would be more realistic— We use exponential fits to 
allow a semi-quantitative comparison of an overlap-only case (Fig. UJ with a situation where 
both effects are present (Fig. EJ)- By considering the gap behavior, we obtain A\ ph 300 
meV and Ai ~ 3 nm for the case considered in Fig. El and A 2 ~ -500 meV and A2 ~ 2 
nm for the case in Fig. [71 Fig. UJ clearly shows the miniband broadening effect as the inter- 
dot distance is reduced, since the states here represented are, strictly speaking, the | el > 
miniband minimum and the \hl > miniband maximum. From the analysis of these results 
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FIG. 7: Single-particle electron and hole energies (ei el> and e| /ll> respectively) and QD gap (ei el> — 
e \hi>) as functions of the distance between dots along [001], in the case of artificially strain- 
unaffected QD's. The solid curves are phenomenological exponential fits, while the horizontal 
dashed lines are their asymptotes. Observe the different vertical scale from Fig. |f)J 

we come to the following conclusions regarding dots separated by distances which are at 
least a factor of 2 larger than the corresponding dot dimension: 

1. The range of the strain field interaction between dots is larger than that of the wave 
function overlap region. 

2. In all general trends shown here, strain effects override direct wave function overlap 
effects, leading to the opposite behavior of the calculated energy variations versus 
distance. The net strain contribution to e| e i> and e\hi> would correspond to the 
subtraction of the data given in Fig. from the corresponding frames in Fig. El Fig. EH 
shows that the electron level downshifts when the inter-dot distance increases, while 
the hole level rises. This behavior comes from the better relaxation of the bond lengths 
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with the thicker GaAs region, shown in Fig. EJ This gives rise to a smaller bond- length 
component of the strain in the ETB Hamiltonian, and then, according to Fig. HJ to a 
smaller electron energy and a larger hole energy. 

3. The miniband width for both the \el > and \hl > minibands is less than a few 
meV, consistent with the largest island spacings considered by Pryor,— where a rather 
different dot geometry and range of distances have been investigated. 

IV. CONCLUSIONS 

We have generalized a previous ETB 2nd nearest neighbors parametrization by Boykin 27 
to include the lattice distortion into the Hamiltonian. We introduced a scaling law of the 
hopping Hamiltonian matrix elements with exponent n = 3.40. We were able to reproduce 
the volume deformation potentials corresponding to the direct (T — T) and the indirect 
(r — L) band gaps for both InAs and GaAs . We have used this approach to calculate 
the electronic structure of a square based pyramidal quantum dot. The comparison with 
previous empirical pseudopotential calculations shows that the ETB model provides accurate 
results for bound state energies and corresponding wave functions. 

The influence of strain on the bound state energies is analyzed. For single dots we found 
the strain increases the main gap by about 25 %. Strain causes the electron states to become 
shallower and the hole states to become deeper. The spatial orientation of the first p-state 
(|e2 >) depends on the alternating interface structures of the four {101} facets, while the 
spatial elongations of the ground electron (|el >) and hole state (\hl >) depend on the 
mesoscopic C^v symmetry of the strain field. 

We have quantitatively discussed the influence of the dot-dot interaction on the bound 
states due to both strain field and wave function overlap by decoupling these two effects. For 
well separated dots, we have shown that the strain field dominates the level shifts, leading 
to opposite trends as for pure wave function overlap, although at distances less than twice 
the dot diameter the latter becomes noticeable. The QD gap between the electron and hole 
states decreases as the GaAs region between dots gets thicker because this allows the bond 
lengths to further relax. 



19 



Acknowledgments 



This work was partially supported by the Brazilian agencies CNPq, FAPERJ, Instituto do 
Milenio de Nanociencias-MCT and by the Deutsche Forschungsgemeinschaft, SFB 296. We 
acknowledge E. Penev for contributing DFT results and A. S. Martins for useful discussions. 



* Corresponding author. Tel.: +55-21-25627336; Fax: +55-21-25627368; Email address: 
rsantop@if.ufrj.br. 

1 D. Bimberg, M. Grundmann, and N. N. Ledentsov, Quantum Dot Hetero structures (Wiley, New 
York, 1999). 

2 S. Ruvimov, P. Werner, K. Scheerschmidt, U. Gosele, J. Heydenreich, U. Richter, N. N. 
Ledentsov, M. Grundmann, D. Bimberg, V. M. Ustinov, et al., Phys. Rev. B 51, 14766 (1995). 

3 Y. Nabetani, T. Ishikawa, S. Noda, and A. Sasaki, J. Appl. Phys. 76, 347 (1994). 

4 J. M. Moison, F. Houzay, F. Barthe, L. Leprince, E. Andre, and O. Vatel, Appl. Phys. Lett. 
64, 196 (1994). 

5 G. S. Solomon, M. C. Larson, and J. S. Harris, Appl. Phys. Lett. 69, 1897 (1996). 

6 H. Lee, R. Lowe- Webb, W. Yang, and P. Sercel, Appl. Phys. Lett. 72, 812 (1998). 

7 J. Shumway, A. J. Williamson, A. Zunger, A. Passeo, M. De Giorgi, R. Cingolani, M. Catalano, 
and P. Crozier, Phys. Rev. B 64, 125302 (2001). 

8 P. W. Fry, I. E. Itskevich, D. J. Mowbray, M. S. Skolnick, J. J. Finley, J. A. Barker, E. P. 
O'Reilly, L. R. Wilson, I. A. Larkin, P. A. Maksym, et al., Phys. Rev. Lett. 84, 733 (2000). 

9 J.-Y. Marzin and G. Bastard, Solid State Commun. 92, 437 (1994). 

M. Grundmann, O. Stier, and D. Bimberg, Phys. Rev. B 52, 11969 (1995). 

1 A. Vasanelli, M. De Giorgi, R. Ferreira, R. Cingolani, and G. Bastard, Physica E 11, 41 (2001). 

2 C. Pryor, Phys. Rev. B 57, 7190 (1998). 

3 O. Stier, M. Grundmann, and D. Bimberg, Phys. Rev. B 59, 5688 (1999). 

4 L.-W. Wang, J. Kim, and A. Zunger, Phys. Rev. B 59, 5678 (1999). 

5 T. Saito and Y. Arakawa, Phys. Stat. Sol. C 0, 1169 (2003). 

6 T. Saito and Y. Arakawa, Physica E 15, 169 (2002). 

7 T. Saito, J. N. Schulmann, and Y. Arakawa, Phys. Rev. B 57, 13016 (1998). 



20 



18 G. Allan, Y. M. Niquet, and C. Delerue, Appl. Phys. Lett. 77, 639 (2000). 

19 S. Lee, L. Jonsson, J. W. Wilkins, G. W. Bryant, and G. Klimeck, Phys. Rev. B 63, 195318 
(2001). 

20 S. Lee, J. Kim, L. Jonsson, J. W. Wilkins, G. W. Bryant, and G. Klimeck, Phys. Rev. B 66, 
235307 (2002). 

21 G. W. Bryant, Physica B 314, 15 (2002). 

22 G. Klimeck, F. Oyafuso, R. C. Bowen, T. B. Boykin, T. A. Cwik, E. Huang, and E. S. Vinyard, 
Superlatt. Microstruct. 31, 171 (2002). 

23 A. D. Carlo, Semicond. Sci. Technol. 18, Rl (2003). 

24 A. J. Williamson and A. Zunger, Phys. Rev. B 59, 15819 (1999). 

25 Y. Fu, Q. X. Zhao, F. Ferdos, M. Sadeghi, S. M. Wang, and A. Larsson, Superlatt. Microstruct. 
30, 205 (2001). 

26 C. Pryor, Phys. Rev. Lett. 80, 3579 (1998). 

27 T. B. Boykin, Phys. Rev. B 56, 9613 (1997). 

28 J.-M. Jancu, R. Scholz, F. Beltram, and F. Bassani, Phys. Rev. B 57, 6493 (1998). 

29 J. C. Slater and G. F. Koster, Phys. Rev. 94, 1498 (1954). 

30 Q. K. K. Liu, N. Moll, M. Schemer, and E. Pehlke, Phys. Rev. B 60, 17008 (1999). 

31 N. Moll, M. Schemer, and E. Pehlke, Phys. Rev. B 58, 4566 (1998). 

32 P. Keating, Phys. Rev. 145, 637 (1966). 

33 R. Martin, Phys. Rev. B 1, 4005 (1970). 

34 A. J. Williamson, L.-W. Wang, and A. Zunger, Phys. Rev. B 62, 12963 (2000). 

35 C. Pryor, J. Kim, L. W. Wang, A. J. Williamson, and A. Zunger, J. Appl. Phys. 83, 2548 
(1998). 

36 Landolt-Bornstein, New Series III/22a (Springer Verlag, Berlin, 1987). 

37 J. G. Menchero, B. Koiller, and R. Capaz, Phys. Rev. Lett. 83, 2034 (1999). 

38 J. G. Menchero, T. G. Dargam, and B. Koiller, Phys. Rev. B 61, 13021 (2000). 

39 P. Vogl, H. P. Hjalmarson, and J. D. Dow, J. Phys. Chem. Solids 44, 365 (1983). 

40 I. Vurgaftman, J. R. Meyer, and L. R. Ram-Mohan, J. Appl. Phys. 89, 5815 (2001). 

41 T. B. Boykin, G. Klimeck, R. C. Bowen, and F. Oyafuso, Phys. Rev. B 66, 125207 (2002). 

42 W. A. Harrison, Electronic structure and properties of solids (Dover, 1989). 

43 P. Pfeffer, I. Gorczyca, and W. Zawadzki, Solid State Comm. 51, 179 (1984). 

21 



C. G. V. de Walle, Phys. Rev. B 39, 1871 (1989). 

S.-H. Wei and A. Zunger, Phys. Rev. B 60, 5404 (1999). 

A. J. Williamson and A. Zunger, Phys. Rev. B 58, 6724 (1998). 

R. B. Capaz, G. C. de Araujo, B. Koiller, and J. P. von der Weid, J. Appl. Phys. 74, 5531 
(1993). 

M. Fuchs and M. Schemer, Comput. Phys. Commun. 116, 67 (1999). 

The three center integrals involve hopping matrix elements between orbitals centered at 2nd 
neighbor sites (A and B) whose values depend on the position of the common nearest neighbor 
site (C) in the zincblende structure. For distorted configurations, the relative displacement 
between A and B has been taken into account in the direction cosines^, while the displacement 
of C with respect to A and B has been taken into account by interpolating linearly from the 
perfect structure matrix elements^, according to the deviation in the position of C with respect 
to the AB line. 

C. Tserbak, H. M. Polatoglou, and G. Theodorou, Phys. Rev. B 47, 7104 (1993). 
L.-W. Wang and A. Zunger, J. Chem. Phys. 100, 2394 (1994). 



22 



